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PREFACE 


A modern high speed digital computer uses many vacuum tubes in 
bistable circuits. These circuits and hence the computer, are susep- 
tible to the failures of the vacuum tubes. Although considerable 
progress has been made towards improving the reliability of vacuum 
tubes, they do not approach that of the later-developed "solid- 
State" devices. In addition to greater reliability offered by these 
devices, the small size and low power consumption is particularly 
enticing where large numbers of  bistable circuits are needed. In 
particular, the development of the ferrites and titanates for high 
permeability and dielectric constants introduced non-linear elements 
in a component size that surpasses the hignly lauded transistor. As 
with other new circuit elements much is left to be done, especially 
in improvement of the materials involved. However, the first step 
was taken by C. F. Spitzer /7/ in the analysis of the two-branch, 
bistable, series resonant circuit, using a non-linear inductance. 
The next move is the analysis of the non-linear capacitor circuit 
which this paper purports to do. The next obvious step would be the 
inclusion of both devices in a single circuit, and a short develop- 
ment of this appears in hopes of beginning the solution to this prob- 
Jem. 

The analysis is to be done with the thought in mind of obtaining 


a set of design equations for a circuit which might be used as a bi- 


zr 





stable network. Both the two-branch, series resonant circuit, which 
consists of series elements in a parallel network, and the two-mesh, 
parallel resonant circuit, which consists of parallel elements in a 
series network, are analyzed. Consideration of triggering Schemes, 
although of great importance, occupies a small portion of this paper 
since it was not the primary purpose to design an actual circuit but 
rather obtain the criteria which might lead to a realizable circuit. 
The production of an actual operating circuit awaits the development 
of materials with lower losses and slightly greater non-linearities. 
The necessary changes and the progress being made on this phase of the 
elements are entirely beyond the scope of this paper. 

The author wishes to express his deep gratitude to Charles A. 
Rosen and Charles F. Spitzer of the General Electric Laboratories for 
their extensive suggestions and guidance in the completion of this 
paper. Without the many "skull" sessions with these gentlemen, this 


paper probably would not have been written. 
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CHAPIER ፲ 


INTRODUCTION 


The circuits to be analyzed in this paper are the two-branch, 
series resonant circuit shown in figure I-la; and the two-mesh, para- 
llel resonant circuit shown in figure I-lb. Both of these circuits 
utilize only one non-linear element. In chapter IV a brief develop- 
ment appears of the two-branch, series resonant circuit using two 
non-linear elements as shown in figure I-2. 

Inasmuch as these circuits contain non-linear elements, methods 
other than linear network analysis must be employed. There are several 
approaches for analyzing such circuits based upon developments by 
Duffing, Poincaré , Kryloff and Bogoliuboff and others /T, 2, 3, h, 
and 8/. The choice of one approach over another depends on the accur- 
acy desired, the accuracy obtainable, and the exact problem confront- 
ing the analyst. The method of Dryloff and Bogoliuboff is certainly 
more general than that of Poincaré; and both are more accurate than 
that of Duffing. On the other hand, the method of Duffing is far 
easier than the others and yet yields the desired information to a 
sufficient degree of accuracy. All of these methods mentioned are 
analytical and suffer from one inherent disadvantage: they produce 
solutions for only the case under consideration without any thought 
as to the totality of solutions available. A qualitative method, 


known as the topological approach, gives an overall picture of the 





ከ 


s 
43 L L | 


E. SIN wt 


QUC 


JQ 


o T 
| 


C= £(9) 


<— 0> 
X 


Figure I-la Two Branch, Loss-Free Series Resonant Circuit 
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Figure I-lb Two-Mesh, Loss-Free Parallel Resonant Circuit 
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totality of solutions but yields no specific information on any partic- 
ular case / 1, 3, and 8 /. 

The non-linear capacitors available at this writing had relatively 
large losses, were dependent on temperature, and had small non-lineari- 
ties. Furthermore, there was no available test equipment for measur- 
ing the component parameters at the higher frequency of circuit opera- 
tion. Finally, the determination of the component parameters was based 
on some rather sweeping assumptions and hence accuracy was sacrificed 
ab initio. 

For these reasons, and due to its simplicity, a modified form of 
Duffing's method was chosen. In order to clearly understand the devel- 
opments of chapters II through IV, a brief presentation of this method 
and the reasoning leading to it is given in what follows. 

Consider the motion of a simple pendulum having no losses and not 
under the influence of an external force. The differential equation 


of this system, as shown in figure 1-3 5 


K+ Jam = O. 16:1 


The solution is begun by using the substitutions 
dx ; ቋ"ሯ . qv d a AY 
BEY] HS VE 


Equation 1-1 then becomes, after separating the variables: 


vdy = -$ anx dx. 1-0 





Integrating this expression gives: 
ve 24 Wax +A, gs 

where h is the energy constant of the system. From this equation a 
plot of v versus x is made as shown in figure I-h. From equation I-3 
and the plots of figure I-l, the following observations can be made: 

l) h >ሪ or V? becomes negative, which defines an unreal- 
istic case: 

2) it Zeche ፦ ; the curves are closed, encircle the points 
v-0, x = 2nTr , n being an integer, and hence the pendulum oscillates 
about the angle x - 2n ; and 

3) if n>34 , V is never zero, and thus the pendulum con- 
tinues to rotate about its pivot. 

The heavy plot corresponds to the case ካ= # , and it can be 
shown that this case cannot exist. However, this value of h divides 
the totality of solutions: those that are periodic and those that are 
aperiodic. 

This procedure just outlined represents the crux of the topologi- 
cal approach. From it there has been determined, once and for all, 
the totality of solutions in the system, but nothing has been learned 
as to the values associated with any rarticular case. 

It seems evident that if a driving force were included, then 
whether or not the solutions were periodic or aperiodic would depend 
on the amount of energy this force adds to the original system. 


Obviously, if aperiodic solutions are allowed in the electrical 
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Figure 1-3 Simple Pendulum 
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Figure I-À Velocity versus displacement curves 
for simple pendulum. 





circuit cases, then, at some time, the currents and/or charges will be 
of sufficient magnitude to destroy the circuit elements. Thus analysis 
can immediately refrain from these solutions and concentrate on those 
in which periodicity is mandatory. 

In order to obtain these periodic solutions we now consider 
Duffing's method which was first applied to the equation bearing his 


name: 


ጋሪ ጋጋሪ le, 0 ee F eut, I-h 
Note that the function of x appearing in this equation resembles the 
first two terms of the expansion of the sin x function. Hence, it 
might be expected that the solutions are much the same as those of 
equation I-l. Such is found to be the case within the region of 
periodic solutions / 1 and 8 /. 
First, assume no damping, i.e. C = O, and write the equation in 
the following form: 
Y uU a NOT X y F ln wt, ፲--5 
The iteration now starts with the first approximation (the approxima- 
tion x =- O being considered trivial), Xu” Aw wt, and inserts this 


into the right hand side of equation I-5. This obtains: 


u 3 
x,» -(&À «3a, AT F)toout - 14 À Coa 300, 1-6 
Integrating this expression twice, there results 
S 
3 As A 
% = zs GA 3ፌዴ4-ፆ)‹= ሬቶ geo Cag UT, ET 


the constants of integration having vanished in order to preserve 





periodicity. 

At this point Duffing made a rather far reaching choice. In lieu 
of the usual procedure of prescribingwand then determining the ampli- 
tude, A; of the fundamental, he prescribed A, = A and thus obtained 

wi= 0, + 20,4 - ተ 1-8 

It is noted that W becomes a Single-valued function of A but that 
A could be a triple-valued function of w . It was for this reason 
that A was prescribed, and hence, w determined. 

The next step of the iteration procedure would be to substitute 
the solution 

x = ACea wt 4 ረ; Ceo. 3022 
36wm* 
into the right hand side of equation I-5 and again integrate twice 
letting the constants of integration vanish. This, then, is Duffing's 
method. For rapid convergence a, and A should be small. 

The modified procedure used in this paper removes the restric- 

tions on 84 and A, but implies that the driving force is small of 


order a, i.e. we set F = a:F. 


3 
This begins by adding ux to both sides of I-5 and obtaining: 
Rue WAL Ya Bist, 1-9 
Again, iteration starts with the approximation Xo= ÁauowZ and substitutes 

this into the right hand side of equation I-9. 


This gives: 3 
X +w z L(w.a,)A-3 34 HAN Ceo w É 


I-10 
- 43 45 Cea 3, 





የሩኪሩቪጂረ፡፡. 





Figure I-5a A, Relationship for Duffing's 
Equation 
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Figure I-5b A Response Curves for 
Duffing's Equation 





Now to perform the next step in the iteration it is noted that the 
coefficient of Cut À 0 or else secular terms of the type Bt Cea wt 
arise. 


This procedure gives: 


(ህ*-ሷ() ጳ - 3 25 4 +፡55 = ዕ 


or à 3 2 A. E 
QJ = O) ት 4 da À ነች E p 
which is precisely the same relation resulting before if 5375 is re- 


placed by F. After setting the awt term to zero, there results: 

ž +wx = -lA m But 
which gives as its solution 

% = A laut + eR. A Cea swt. 1-12 
As in Duffing's method the choice is to make A, = A. The procedure is 
then to substitute I-12 into the right hand side of equation I-9 and 
repeat the steps above. The iteration is continued until the desired 
degree of accuracy is obtained. 

It is instructive, and will be useful later, to plot equation 

1-8. Although w* is the dependent variable in this case, it is plot- 
ted on the abcissa in order to show the resemblance of the response 
curves to the linear case. The plots in figure I-5 are for the loss- 
less case, while those in figure I-6 are for the case where losses are 
considered / 8 /. It might be added parenthetically at this point 
that the plots in these figures bend to the right or left according 
to whether a, is greater or less than zero respecitvely. If a, is 


zero, the linear case results and the curves are upright. 
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Figure 1-7 Q-E Relationship for A Non-Linear 
Capacitor 





Figure 1-8 Nest of Q-E Curves for Varying 
Field Amplitudes 
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If a single series circuit were to be driven with a constant 
frequency source whose amplitude could be made to vary, then operation 
would take place along a vertical line such as that show in figure I-Sa. 
It should be noted that as the amplitude of the driving force is slowly 
increased a point is reached where there are two possible values of the 
response, A. since A represents charge in the electrical case, and 
current is the time rate of change of charge, when this critical 
force is reached, the current is said to jump. If, instead of a 
single series circuit, two such circuits are placed in parallel as in 
figure I-la, and a constant current source used, it would seem plaus- 
ible that when the critical force level is reached, one branch current 
would "jump" up while the other branch current would "jump" down. Such 
an event does occur. 

The question now arises as to what form f(q), shown in figure I-l, 
shall take. If a charge versus voltage plot is taken on one of these 
non-linear capacitors, a curve such as that shown in figure I-7 results. 
This particular curve is for one value of driving voltage. In figure 
I-8 a set of such loops is shown for various driving amplitudes. The 
heavy line connects the loop tips. Obviously, no simple relation de- 
fines the Q-E curves of either figure I-7 or 1-8. 

oince the real interest is in the value of charge which can cause 
bistability, it would seem evident that the large signal loops are of 
more importance. This is due to the fact that the charges must be 
large enough to cause an excursion into the non-linear region at the 


time of "jumping" into bistability. 
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If a single-valued function is desired, it is assumed that the Q-E 
relationship follows the mid-points rather than the upper or lower por- 
tions of the loops. This tends to average the characteristics and has 
the advantage that voltage is zero when the charge is zero. 

The assumptions made this far are the primary Jimitation on the 
accuracy of any solution. It would seem rather irrational to 
choose an involved functional relationship however accurate it may be. 


Thus the defining relation is chosen as: 
ሷ 3 
— + 5፡7 3 1-13 


an odd function which is required for the midpoint loci. 
The necessary information has now been obtained in order to 


commence the analysis of the bistability problem. 
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CHAPTER II 


TWO- BRANCH, SERIES RESONANT CIRCUIT 


It is the plan of this chapter to investigate the two-branch, 
series resonant circuit. Investigation first concerns itself with 
the loss-free circuit in order to determine the basic relationships 
for bistability. After these conditions have been studied, emphasis 
will shift to the low-loss circuit, in an effort to determine the 


effect of resistance on the bistability relationships. 
SECTION I: THE LOSS-FREE CIRCUIT 


The equations of the circuit of figure I-la can be written as 


follows: 


e = LE ት ር + = E 1 


ay, ያጋ 
ር me = €, js A qu 


Subtracting equation II-lb from II-la, there is obtained one 


of the basie equations for the circuit: 


L Dh) GS 


Lh 


E 
E 
Po 





The other basic equation for the circuit is: 


bores = 


As in chapter I the defining relations for C, and C, are 


chosen as 
Es = A: 1 ት Oa E an 
C, = Ay fa + 43 ዴ : Il-hb 


oince it was decided that the circuit should have a constant 
current source, the impedance seen by the generator voltage, Es 5 
must be that offered by (4, otherwise the varying inpedance offered by 
the two non-linear branches would vitiate the constant current hypoth- 
esis. Thus ig leads eg and hence, since it was tacitly assumed that 
there were no losses, the first iteration begins with: 
dio = ዳ tol II-5a 
= Ban we, II-5b 
where the second subscript denotes the iteration by which the value 


was obtained. 


Equation II-2 is rewritten as follows: 


( .- =) +ሠሂፃሜ) == ኣሚት”) —L (Cem EEG 


The next step is to substitute equations II-5a and II-5b into 
the right hand side of equation II-6. This results in (see appendix 


III for expansions of powers of trigonometric functions): 
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#፦ ) = ሠ'ቆይ-ይ) = [wr e. )(A-B) - 24 Bonot 
ት ee B uengut 11-7 
The conditions for periodicity require that: 


( በ ፦ ቷ ነረሉ- B) - 343 (43- B?) — 2) TIS 
P 2L 


Disregarding the trivial case of A - B and using equations 


II-3 and II-5 there results two equations relating A and B: 


Y=o= —-6 + A +65 m 


Y=z0= -5 + m(A"+AB TB. 11-10 


where 


UD) = 
S= wi | mu 
Q. 
Wo = T ) Rey? 
1 ZQ i 
m = >, II-13 
ቃመ * 


ሪ is a straight line whereas 4 is an ellipse. Using the rota- 


tion matrix 


A ሽ = - 72 A TU 
IL 
D ez IIA 
equation II-1LO becomes: 


ይጆዕ = A» + I” +67 B x 


y 
lo 


55... 





The half axes are 

Us = y 33 ) and II-16 
Zs 
N 


2 II-17 


For the two-branch series resonant circuit the second iteration 
must be performed in order to obtain bistability conditions. The only 
information obtainable from the above relations, is that in order to 
have real solutions of Ya ) S > O and hence Ww E We, where 
y means: "must be greater than". This condition can be deduced 
intuitively as follows. For bistability A ae B and hence one of the 
two identically designed branches must offer less opposition to charge 
movement than the other, due to the non-linearity of its branch capaci- 
tance. This requires one branch to be "tuned" nearer to U than the 
other. Now for a capacitor defined by an equation such as II-l, it 
is found that the effective capacitance presented to the fundamental 
voltage decreases as the amplitude of the impressed, fundamental 
voltage increases. Hence as the amplitude of voltage increases, thc 
resonant frequency of the branch must increase. Hence the circuits 
must originally be tuned below the source frequency i.e. WI>Wo- 
Thus, what was intuitively believed has been proven by the inequality, 
EC 0. 

Solving equation II-7, after the condition of periodicity has 


been invoked, there is thus obtained 


rie (AB) wi - As (Baila, 8 


እ | 








The choice as explained in chapter I, is to make A =A and B, =B. 
Equating the portions of II-7 identifiable with q} and q respectively, 


there results 


fu = Acin wt - DA” ac 3wZ 11-194 
Qi = R aum wt DEEE ፲፲-195 
where 
as M u 
32102), ጋ#ሃ/57ፇ/3 " 


Substituting these relationships into the right hand side of 
equation II-6 there results (see Appendix II for details of the 


DO: 
(Q, 8), - (2-9) - 
[6-2 4-B) - 30s (4*. 85) 
“23 (ep) da Eg) anet 
+ (RPG 8) + Soir ዕን SEDE on, 302 
- ETE ER ን” 
+| Boa, Az 
+ N (ALE ar Jwt., wa 
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Invoking the condition of periodicity and disregarding the trivial 


solution, A - B, there results: 


2: A ሃ 
ህ3፳ዕ = - 8 + ሰሓ (ዳኙ 48 ተፀደ*) + nn D(A + ÆB+AB "+48 +B 2 
+ 2 mn D'CACHASB+ AB BT MB GAB +B), 55 op 

Using the rotation matrix of equation II-1), equation II-22 


becomes 


"t 


0=-2 +5 (3% ?+8°) + D (So "+ +ያ) 


R 
+2 ( V 6435 ሪ +2/=2'፦ሪ*) MEOS 


Va 


This form is still not too instructive. However, another trans- 
formation reduces the form to where the desired results can be obtained. 


Letting ez Cie à and A =P amó , II-23 becomes: 


| soc A +4 fop* +246" 246 Mí cna 


where £@) - 2+ C020, II-25a 
99) = 5 ፉሪርሪሪሪፓፀ -ሠፓርፎ ፉሪ II-25b 

and h(6)= m 04006) 4- 18 CW 26 - 40°26] 
2E 7 ዕ ሖኃ 8 : ፲ ፲-ፎ5ር 


For any value of O, % . has only two finite, real roots (posi- 
tive and negative values of the same magnitude) and hence represents 
a single closed loop.V/, Y and Ye are shown in figure II-l. The 
developments of these curves and the proof of the single root magni- 
tudes are shown in Appendix II. 


From figure II-1, the mode of operation may be explained. For 
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small values of Ez, and hence Qg , the two series branches behave as 
linear networks and thus A = B. As Es is increased, Y moves out fr 
the origin parallel to the axis. The operating point remains on the 
intersection of # and the line, A= B. In Appendix II it is shown 
that for O - O degrees, the slope of the tangent of e is parallel 

to the 2 axis. Hence X will be tangent to E at the operating point 
P. When Y becomes tangent to Vo ; there are now, two other operating 
points available, Po and P3- Although no mathematical proof is given 
here, it can be deduced from the work of Stoker / 8 / that the area 


around P, is unstable. Operation, therefore, "jumps" to either P, or 


2 
P3- Ihe reasons leading to this deduction of P, being unstable are 
given in chapter III for the first iteration of the two-mesh parallel 
circuit. The reason for postponing this discussion until then is 
apparent from the discussion itself. 

A circuit in which the controlled switching of the operating 
point from P, to En: vica versa, is obtainable constitutes a trigger 
pair. 

The reader might now infer that the firsi iteration will not 
predict bistability for the two-branch series resonant circuit.  How- 
ever, in chapter III it will be shown that this is an unwarranted con- 
clusion; the same value of E, which will cause bistability will be 
predicted by two different approaches. 

For operating point P} in figure II-l, the following values, 


developed in appendix II, are obtained (9 - 00): 


2l 





y 
= 258 | A * 


; da / o 
AU = = #9 S 4 22. I1-27 


S z/ 
Ds 


Ls 


gv S ya e 11-98 
2o? 


ዘ 


UJ Os II-29 


Obviously, if Y moves out from the origin so far that is fails 
to intersect 4 , bistability fails to exist and both branches revert 
to equal currents. Form figure II-l this appears to occur when is 
tangent to fs at about OG - 209. For this condition, the value of He 
becomes 

Is= w = we, Les 20° 4,010 20°] 
= #2 wp aon 65° 11-30) 

Greater accuracy can be obtained by finding the value at which 
ds vanishes, but this hardly seems like a reasonable procedure in 
view of the original assumptions. 

It should be added at this point that after bistability occurs, 


there is no assurance that increasing Eg will increase I This is due 


s 
to the phase relationships possible between À and B. If the two opera- 
ting points are located in the first quadrant of figure II-l, then A 


and B are in phase. If the two operating points are located in the 


second or fourth quadrants, then A and B are 180 degrees out of phase. 
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SECTION 2: THE LOW-LOSS CIRCUIT 


Attention is now given to the case where the two branches contain 
losses as shown in figure II-2. Equation II-6 is rewritten to include 


the effect of resistance as follows: 


(9-92) + e 70) = W(9,-92) -£ ( 7 -9.) = Ce 0). EY 
Assuming that R is small compared to Xc} > then I, will still be 
primarily determined by the impedance of C}. For the remainder of the 
discussion in this chapter, it will be assumed that AN» B. Hence, I, 
is very nearly in phase with I,. Since Is will lead Eg by 20 degrees, 
the first iteration begins with the approximations: 
dr À Lin wt o 
$2.7 Barn (eet -4) 11-32 


Substituting these values into the right hand side of equation 


II-31, there is obtained thereby: 
Gr) + weg.) = Lur- EMA- 885] ac ue 
- [(w* &)B - 32 B. daz. (wt- p) 


LU » 
= EW Acsawt + £4 B tse lwt- y) 


ı 93 A; amwt - Qs Bun 3(wt-p). 
Eu re # p) 11-33 
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Figure ፲፲-ረ 


Two-Branch Low-Loss 


Series Resonant Circuit 





Figure II-3 


Phase Relationships for Two-Branch, Low-Loss Series 


Resonant Circuit. 
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In order to preserve periodicity all of the fundamental component 


terms of W must vanish identically. Hence there results: 


| -SB +mB* lan (wt-y) 4 wk B t (t -p) 
, [E sA tma ainw ዣ ፈራ AcsowiE, 11-31 


Expanding the terms involving (ሪ/ድ-ዖ) and equating the coefficients 


of sin wt and cos wt respectively, there results: 
E gs ሠ LAY = |- 3 R 
EsbtmBlwop + YEBeiny= ESA +m], 11-352 


= US 
WE B cea p- EsB+m8?] any = A 11-35» 
The phase diagram appears in figure 11-3. Since R was presumed 
small, y is very nearly zero degrees and hence cos ሥ mw le 6 “Usang 


this fact and dividing equation II-35a by II-35b, there results: 


LA B den p — ESA+mA - ሁ-5 ጸ+/ ፳5] 
=|-sg+mB any WRA-WÄB 


or 

ሙ 5 + m(A+AB+ B*)- = e es 11-36 
If R is allowed to vanish, the first iteration, loss- free case is 
immediately derivable from the above.  Experimentally ^m B^ is of the 
order 8 x Aone Whereas S is unity or greater, so that the assumption, 
SX/mB, is not too unreasonable. 1f it is allowed that s >> MB“ then 


II-36 reduces to 
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/ a 2 በ 

=0= A +AB+B—2 
X +AB+B — =, a. S, <= 
Now VU, will not be a real illipse, and hence bistability cannot 


result, unless 


as 
À < um ? 11-38 


which places an upper limit on the losses, 


For B20 11-37 gives 


| = _ wi 
== fm T asm ' ገ 


If R is allowed to vanish, the first iteration, loss-free case 
is obtained, as can be seen from equation II-10. 
From figure II-2, it can be seen that: 
eu LI +R 9, ES II-10 
or 


E AS ESA mA Tac ut + WEL coat f, 11211 


The magnitude of E, can be found by first finding its square as: 


2 ጊ 274 Le 2 
= = a [-$AFmA ] + A A X TO 
Substituting the value of A from II-39, and upon simplification, 


there results: 


E = ፊያነሃ/ጹ - ጄድ Y, are”) ee 





N G3 Sm 
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Since = 








II-l 
and, since 
> S E 
Il=60B= &, 
PRS) / = — Qik 
Se QS (3 ea ; 
11-15 


The ratio of A to B (i.e. the ratio of the larger to the smaller 
charge magnitudes) is defined as the switching ratio,4A , and can be 


obtained from equations 11-39 and II-15: 


ሷ,= 
Ac 7t | 11-16 





Now from equations II-la and II-lb, remembering that B470, there 
results: 


E = AAr A, ace 


ሁኒ; = 4,6. II-lTb 
Since the usefulness of this circuit is predicated upon a suffic- 
iently large difference between the two branch output voltages, it may 


be desirable to solve for this voltage difference: Hence, by sub- 


2T 








traction, and by using equations 11-39 and 11-15 there results: 


2 pe 
E Es lo GRA PE, a 


which for the loss-free case gives 
ሃሯ 
a ee 0, (1+S) ሥታጋን e 11-19 


Thus it can be seen that although the switching ratio, À , is 
independent of the degree of non-linearity, m, that the voltage diff- 
erence 38 inversely proportional to its square root. It should be 
noted, however, that À is not the ratio of E to BJ, for the circuit 


under consideration. This value, still independent of m, is given by: 


E, . ASI 5 — WER 
E ^ er [ITS ays 


z ሷሳ [|+5 - Ze 2 II-50 


The smaller the value of m is, the larger the potential difference 
between the non-linear capacitors is for a given A . This follows 
from equation II-l9 and the fact that À is related to the parameter, 

s, by equation II-l6. This potential difference, however, is obtain- 
able only by making A larger than B. put since the switching ratio, 

A > has been prescribed, both A and B increase, with B remaining in 

the linear region and A in the non-linear region of the Q-E relation- 
ship. 

However, although A is independent of m, it must be remembered 
that too small a non-linearity results in currents of such magnitude 
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that they may destroy the components before bistability is reached. 
A11 of the above relations were derived from the first iteration 
of the circuit with small losses. The second iteration is now com- 
menced in order to find what modifications, if any, need to be made 
to the first iteration relations. 
After the periodicity conditions have been invoked, equation 


II-33 becomes 


G si 9) 22 Ag) = 2: Tun ያ2ሬ/ሯ 


ቺ m Eo 3(at- p), 


TI-51 


whose solution is 
| = 
9-2) = (A-8)atwzZ - DA am Ft 
VS, HET - 10) y I 
Again choosing ÁjzÀ and BizB there results, for the second approxima- 


tions: 


f" = Agen ሪሪሯ -- DA Sco II-53a 


f = Bamwr === BEP 3(WwZ- y), II-53b 
Substituting these expressions into the right hand side of 11-31, 


gives: "ee 


Qu je), + ሚጎ“) = 
WA A-R- ZER) 3 PAE jam? 
፡፡ መመር wt + S s DB n jut 3/0) 
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— WE (A-B) ዕዕሬ ሠራ 


Ls 
+ «ና € 4 a 7 ፈ 9 , ሪዐ “¿Mu 7 » ) II-5l 
where the dots refer to terms involving 3WL, 548 UT end Jwt, 


The conditions for periodicity become: 


—UR (A— m E 
A (A B) + $43 DB 43-0 11-55 


[(Lo^- %)(A-B)- 3as(42 82)- 59 ዕፆኝ ድፍረ ) 





_ 343 n2/4’_EN] = O. II-56 
343 n'(4"- B^)] 
From II-55, the phase relationship is determined as 
| - ሄሪያ (4-8) 11-57 
Thus, since rr B Y = | 
RS a, .B' II-58 
— Ya (A-B 
Solving for A-B results in: 
= 2. 303 BT 
4-5 = “መሀሥራ 
.| < 3% 8’ 
e QS mB' 
A= irs ተ! 
፲ 2 


Assuming that equation II-ló6 does in fact give a reasonable 


approximation for À , there results, upon substituting into II-59 


ara ው 
= S-mMB65' + 11-60 
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Here again since bistability implies a value of À greater than 
unity, a small value of m may necessitate a value of B so large that 
the components themselves might be destroyed. 

Substituting equation II-57 in equating II-56 and by using the 
first two terms of the expansion of Cea 3 D; there results, assuming 


to be small: 
S (A-B) — ^m (A*- B*) - mb (A B^ Sd Ron 3p-B)) 
— amab*(A-B?)z 
II-61 
Collecting terms, excluding the trivial solution, A=B; applying the 
rotating matrix, II-1l; and using the polar coordinate form, there 


results for O - O9: 


ጄ = Ze 1168 andv - #7"? 11-62 
&o° 


/ ጋ 05 ¢ 
= TT M 7 ጅ ብራ “ድ E 2034 ADS, a) yy 63 


Except for large values of R and NL , V” is approximately equal 
to dk as given by equation II-27. 
Briefly sumarizing these results: 


1) in order to have bistability 


= EX 750) 6 Paci I > 11-26 


where AL = tdink 2) mits ድመ 7 11-27 


2) the switching ratio is given ጅን 


ኢር : 
A= CGR 2 E 


dl 











3) circuit losses should be such that 


Quas — 
K < ELOD II-39 


l) the difference of non-linear capacitor potentials is given by 


ኑ፦ " = a (4 S)) A ዘ] II-19 


5) the ratio of the two non-linear capacitor voltages is: 


E A (19 S = Zh); and 1250 


6) the non-linearity should. be such that 


S > > 
E 2 À >| 11-60 





does not require B so large that the component ratings are exceeded. 
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CHAPTER III 


TWO-MESH, PARALLEL RESONANT CIRCUIT 


The plan of this chapter is to develop relations for the two- 
mesh, parallel resonant circuit similar to those discussed in chapter 
II for the two-branch series resonant circuit. It will also be shown 
that the requirement of a constant-voltage source in the place of the 
previously used constant-current source gives rise to a different solu- 
ticn curve relating, A,B, and fs, the source voltage. 

As in chapter II, the first consideration will be the loss-free 


circuit and then the low-loss circuit. 
SECTION 1: THE LOSS-FREE CIRCUIT 


The equations for the circuit of figure I-lb can be written as 


follows: 


L = Lede P ap III-la 


==. yf e, de + f: III-1b 


ፈ 





Subtracting equation III-lb from III-la there is thus obtained 


one of the basic equations for the circuit: 
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7 fle-ejdt + LGS) - o. uma 


Since there are no d.c. sources in the circuit, no loss of 
generality occurs if equation III-2a is differentiated with respect 


to time. This differentiation results in 


Gf + P (e -€,) =0. III-2b 


This equation is the same as was obtained in chapter 11 for the 
two branch, series resonant circuit. 


The other basic equation of the circuit is: 


CE Cou (Ce. ንክ 
As with the series circuit, the defining relations for ©, and Ss 


ares 


a Qi 1. + âz 7. III-la 


C= aR + às; ፲፲፲- ከ5 
In lieu of the constant current source, a constant voltage source 
is desired. Hence C, should be fairly large. Thus, iy is a small 
quantity. For lossless circuits, C, and (C, both lag ig by 90 degrees 
if W>Wo. The charge will be in phase with the voltage and the 


first approximations can be written as: 


dio — A Coo we In. 


G20 = Benw wl. III-5b 


Proceeding as in chapter II, the two solution curves relating A 


and B are: 


3h 





02-45 ፦ ጋሚረዳ” ፉ4ሮ > 2”), III-6 


S == = + (A+B) + (AH ES ፲፲፲-? 

It should be noted that the X curve is common to both systems, but 
that the straight line, y , has now become the cubic ፊ . These 
curves are shown in figure III-l, as derived in appendix II. The fact 
that 5 is a cubic as shown, makes it possible to demonstrate bistabil- 
ity from the first iteration. The operating point for low values of 
E. is on the intersection of ps and the line, A-B. ለ8 E, is increased 
fj moves out from the 6 axis, its asymptotic line remaining parallel 
to the A axis. When Yy becomes tangent to A at Pj, the operating point 
"jumps" from P, to P or Pa Mathematical proof that E is an unstable 
operating point is beyond the scope of this thesis. As stated in 
chapter II, however, the instability may be inferred from the work of 
Stoker / 8 /. The reasons leading to this conclusion of instability 
will now be shown. 


In chapter II it was shown that the half axes of the ellipse are 


EL / $3. ) and I11-0 
ሪሪ = 7 m e III-9 


The first of these represents the "distance" from the origin to 


P. Since A-B at this point: 


A = f zmn ° 1II-10 


Substituting this value into equation III-7 results in 


SE 








Es= 20,0 + SEE jan 


As suggested in chapter II there is another method by which 
equation III-10 and hence equation III-11 may be derived. It is the 
criterion of this method which leads to the conclusion that, P, is 
unstable. In chapter I it was shown that the frequency for the loss- 


free Duffing equation was: 


2 
Wut =a, + 2-05 4 -ፓ . I 


Since Duffing's equation is the same as the equation for one mesh 
of the circuit under consideration if F is replaced by the voltage, መቹ 
across the mesh, the frequency is given by: 

wt= a +33 A =. Tee 

The locus of the points where the curve of Fig. I-5b has vertical 
tangents can be obtained by differentiation of equation III-13 with 
respect to A and setting dwd A equal to zero. This procedure results 
in: 

wW- h -$ls A 20 ral) 

The value of A, for a particular W, which lies on the locus 

given by III-Ll;ı is found by solving equations III-13 and III-1h simul- 


taneously. This results in the value of A, 


S 
A = V 357° Tr To 


Stoker has shown that the region bounded by equation III-1) is 
unstable / 8 /. It was assumed that both meshes of the two-mesh cir- 


cuit were identically designed, and that the circuit was supplied from 


ጋ” 





a constant voltage source. Thus, when A reaches the value given by 
III-10, only one of the mesh voltages can "jump" up. The other "jumps" 
down. Hence, Pz is unstable and operation "jumps" to either P5 or Pz. 

The value of E, given by equation III-11 has another interpreta- 
tion which may be obtained by study of figure III-1. If E. were larger 
than the value given by equation III-11, then there is the possibility 
of having two sets of bistable operating points, since A will inter- 
sect ይ in four points. Thus for no more than two possible stable oper- 
ating points 

ሠ 20, (1 +2 )/7ሯጁ ፲ፕፐ-35 

This apparent contradiction between equations III-15 and III-11 is 
explained as follows. Bistability results when Es reaches the value 
given by equation III-ll. However, as the amplitude of Is is increased, 
there is no assurance that the value of Es will increase due to the 
phase relationships possible between A and B. As stated in chapter 
II, operating points lying in the first quadrant represent the in phase 
conditions, whereas those operating points lying in the second or 
fourth quadrants represent the 180 degree phase conditions. However ,if 
Es can be caused to exceed the value given by equation III-12 after 
bistability has been reached, then four stable points may be possible. 

As Eg is further increased, there is a value beyond which bistabil- 
ity will cease and reversion to the A=B condition occurs. This value 
of Es is determined by finding the common point of tangency between % 
and # besides the one already found at Pj . For % and % to osculate 


at a given point, the following relations must be satisfied at the point 
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- As 
ጠ amt III-16 


= 24 III-17 


when Es has been adjusted so that the point of tangency lies on (Z2, 
then both the necessary and sufficient conditions have been met. 
Applying equations III-16 and III-17 to equations III-Ó and 


III-7 there results, respectively: 
| (. 
/ + smA = m 24 +8), III-18 


2 
/ + ምጋ ረፓ = Wr) (28 AA), III-19 
Division of these two equations by each other, and simplification 


gives 


3n (4287 )+ /ሪጋጋያ-///ሪ-2) = 0. III-20 
The solution AzB is trivial for the present development and hence 


there results: 


2! ar N 
W=o = A 34B*8* - E — 
Application of the rotation matrix 
EE 
^ = E M III-22 
B We Wi 
to equation III-21 results in 
ንንን ታሪ ርር Io. 
CR oS ማራ 


a hyperbola with half-axes. 


lza ES a. III-2] 
LE — pu ' III-25 


By 





The ellipse, 4 and the hyperbola, £ will intersect if the half 


axes obey the relations 


// e Ye 25 
ሺ ሠ ሪክ. = 5, m < Sn * III-26 


From this inequality another criterion for bistability results: 


5>ሙ, a 1299 
UJ? >b wi. 111-28 
The two criteria for bistability, i.e. the required values of s 
and Es, thus far derived have been obtained from the first iteration. 
The second iteration should give corrections to these relations. 
Proceeding as in chapter II the two solution curves thus obtained are: 
f3 =0 = "4 H(AAB IE) ANAL ABI B SAB AB 
A 2214 6 4A BHATE ZA BA BUB EU III-29 
GEO = -5È » (A9B) vm (A948) + AGED 
+ 2mD° 4767).  1r1-30 


The curves are shown in figure III-2 Continuing the reasoning of 
the first iteration procedure, bistability occurs when # becomes tan- 
gent to A at P,. Applying the rotation matrix, equation III-22 to 
equation III-29 and equation III-30 gives, respectively: 

A -0-—R +t x47) + lon’ rk v) 
3= 07 m a E # ፖ= 
2. 
= /%ረ* %° 7226 : 
2 (7« +5578 7 2/008 78 ), 111-31 
= E, ድ E Â 2 
= z =L 2202/5, ,A ,2 y 
he =O Fe to +00 (X74 Boye?) 2a spé) 


2 = 
-+ EE 42/2 ፉዱፍናሪያች ይነ +). III-32 
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As before, the solution curves are now brought into the polar 


form through the equations: A Pod and, z^ à . This results in: 


%=0= == "dll yon. Ap A 2%) > 111-33 
EE -85 + 208092 map = 


SnD 5 2 
TE kaye 724) 2” TEE 
where, f(0), g(0), and h(0) are as in chapter II; and 
4 (0) — 3022 — 205238 III-35a 
ሪ /#) — 04 — ftoa!0 III-35b 
LG) = 00/4000 0 -28 C0 2 > ያረኃሪ 22. III-35c 
Figure III-2 shows a plot of these curves. The curve £ differs 
The curve 


but little from AL, except for large values ofan e 


Y, shown is therefore an approximation. 
The value of Es required for bistability is obtained by substitu- 


ting into equation III-3 the values: ÿ= © degrees and 
111-36 


e Pa ን 
Where DIOE El ain (2165 280$) p III-37 


This procedure results in 
E - À A + M7) Z4 mò Sy 20207? 
Serta C ብ 12/7), III-30 


where ሪ is given by equation III-36, 
There is no easily derived correction for the required value of s. 


Application of equations III-16 and III-17 to equations III-29 and 
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III-30 results in an equation which is very tedious to solve.  Con- 
sidering the approximations first made regarding the Q-E relationship 
of the non-linear capacitor, the solution does not seem worthwhile. 


SECTION 2: THE LOW-LOSS CIRCUIT 


As in chapter II, attention is now given to the low-loss circuit 
as shown in figure 111-3, with the phase relationships shown in figure 
III-4. Equation III-2b then is rewritten as: 

ራ5 de - E m s 9 ስ " : 
(6-9) rg 3 > «(4-2 )-1:(6 -&) -1 (e-&) 111-39 
Assuming R to be large compared to XC4 > the first iteration commences 


with the approximation: 
90 = Alaa wt III-h0a 
ho መ ዕሪሬሪ/ ምሥ $) III-NOb 
Substituting these into the right hand side of equation III-39, 
and involving the conditions for periodicity there results: 
Cw- £ A- Zup hemut + Y Ja, A PLA) ont 
- [ር -9)።- $283 ]ece (tut +8) 
- 49 [a B2 6,8] o^ (224) = O0. 1a 
Using the expansions of cos(wi4$) and sin (wÉ-* $) and equating 


coefficients of like functions there results, since cos £% /: 


% (a, B+24 Bain S= (wa YA-B)- 25 (A2B3) ፲፲፲-- ከ25 


t ANB 34s Lans=-2 la 4-8) +30 (4283)]. 11-12 
Dividing equation III-l2a by III-l2b and subsequent simplification 


results in: 


h3 
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Figure III-3: Low-Loss, Two Mesh Parallel Resonant Circuit 





Figure III-l: Phase Relationship for Low-loss, Two-mesh 
Parallel Resonant Circuit. 
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HE Me 4 m (A48 9 8?) 
=i] 2: ኣመጽ ንጋሽ ን ን ርን Tira 


ENTE. Hang? S- m (AFAB +8?) 


2 2 
Since in practice, “mB ` is of the order of 8 x 10 2 it can be 


neglected in comparison with unity. Thus equation III-13 reduces to 


NL us ptam (AAB +E) ፡ III-lji 
Wu? L* " 5- ^n (4*4 AB 757) 


Solving this equation for A ABB? results in: 
2/ > 
2 Bt = $% A cS 
4 #48#*7ጽሮ ፡ፋ፦-፦ መጭ ^ III-li 
wm (Ss _ wid ) 
ጭማ 
Clearly this reduces to the loss-free case as R becomes infinite. 
If the right hand side of equation III~h5 becomes negative, the equa- 
tion fails to define a real ellipse. Thus in order to permit bistabil- 


ity, the following inequality must hold: 


wh, 
pom s III-h6 


For £220 , Ais taken from III-l5 as: 










si. ሟኗ 
E A o = 
mis- ag") III-h7 


Solving equation III-la for igo, and using the fundamental com- 


ponent only results, after simplification, in the expression: 


፲341-ከ6 








—————— — 


oince B is small, the admittance of the associated mesh is: 


we”, ር ው 2 
e | Re = D +s* INEO 


Also, since B is small, there results from equation III-lib: 
lo 


E. > 0,8 = Ya > or III-50a 


luc IEEE LI 
2 2 y Y 
- ጩ (S41 ys ES p^ 
D = EF = a TII-50b 
( 5 ) /፤ሥ/=፦ <= 
Again, defining the mes ratio,Z\ , as the ratio of A to B 


i.e. the ratio of the larger to the smaller charge amplitudes, there 


results: 


A = ie) Sis Wh) 
~ OL ።። | ፲፲ ፲=51. 
Now, from equations III-la and III-lb, and setting BX’O , there 
results: 3 
E = &,À + As A III-52a 
aF Qul. III-520 
Since the usefulness of this circuit is predicated upon a suffici- 
ently large difference between the two mesh voltages, it may be desir- 


able to solve for this voltage difference. Hence, by subtraction and 


use of equations III-l7 and III-50b, there results: 


wè $ 


A, |] Ye = +0 ት III-53 
RIS- Y) ms) 


l6 











For the loss-free case, this reduces to: 
mu ብ Sy <. 
3 | = Q, (I+ 22 y 2= . III-5l 


Here again, although the switching ratio is independent of the 
degree of non-linearity,/», the voltage difference is inversely pro- 
portional to its square root. 


Ihe ratio of the mesh voltages is given by 
E Ss + uy” 

S =A(1+ 4 se | IT1-55 

R 


Since the second iteration of the low-loss circuit produces 





algebraic equations which are not readily solvable and the first, 
iteration of the low-loss circuit has become rather involved, it 
hardly seems reasonable to attempt a second iteration in an effort 
to obtain more information. 

In Summary, then, the following conclusions apply to the two- 
mesh parallel circuit: 

1) The driving frequency must be such that for bistability to be 


possible, ቿን 6 
2 LII-2 
wi > ws; 


2) The minimum value of voltage required for bistability is, 
from the first iteration, 
Sy Tea 
em = 30. (1 2 I/M ) 
3) The losses must be such that 


t -h6 
R > ግጅ ፲፲፲.-ከ 
S J 





hl) The switching ratio is s 
PERTE 
_k (Ss - Se) 
= GE S+] ) 
5) The voltage difference is given by 


ጊ è " n 
Et ET #; ፡ aj- gin 7, 7 “5- 223 vn 


6) The value of E, required for u from the second 





tite ol 


iteration 

-0 = $^ 2447 
Es = & (0 rome + PP "+ 2 ሥሠ ) ፲፲፲ --36 

where y 
254 hé ton « 2417] a ፲፲፲ -36 

2 
© J / IAS 4-805 

and fr = + ind ዶ %>/ PESO 
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CHAPTER IV 


USE OF TWO NON-LINEAR ELEMENTS: 
FIRST ITERATION 


In chapters II and III, the two basic trigger-pair circuits 
utilizing a non-linear capacitor have been briefly studied. The two- 
mesh series resonant circuit, using a non-linear inductance, has been 
analyzed by Spitzer / 7 /. Intuitively it is felt that the inclusion 
of both types of non-linearities would reduce the currents required 
for the same switching ratio. This chapter will investigate this 
possibility by analyzing the loss-free, two-branch series resonant 
circuit using two non-linearities as shown in figure I-2. Due to the 
complexity of obtaining a solution and the inherent inaccuracies pro- 
duced by the assumed current and charge functions, the analysis will 
be confined to the first iteration. 

For reasons given in chapter I, the current-flux and charge- 


voltage relationships, are chosen as: 


SR) mé c4 +4 p7 IV-1 
#(ያ›።* aug +03 7° IV-2 


Furthermore, since á 
= [= dt IV-3 


equation IV-2 can be written as: 


ከ2 








e = afs fó 4 [448] 
ms DET "Af AT IV-l 


One equation for the system is given by 
N ELBA) WAS = o, IV-5 
whereas the other equation is: 
ty EL, = መ IV-6 
The branch currents and voltages are from equations IV-l and 


IV-2 respectively: 


3 
2464 4. Z, IV-7a 
d 


< + 43 D IV-7b 
a, f, i A. IV-8a 


= IV-8b 
C= E+ asd 


Again, the capacitor, C41, is to be the primary impedance offered 


La 


0 
| 


to the driving voltage. Hence Ig leads Eg by 90 degrees. For the loss- 
free circuit, ihe branch currents will be in phase with the source 
current, Is. Thus the branch fluxes will also be in phase with Ig, 


and the first approximation is taken as: 


Pa = A Cod cuz IV-9a 

Pro = 6 Ga e, IV-9b 
Equation IV-5 can be differentiated with respect to time without 
affecting the solutions since there are no d.c. sources. After per- 


forming this differentiation, the equation is solved for (2 = N )s 
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In the resulting equation, 4*(F-%) is added to both sides, and there 
is thus obtained: 
6-4 ow G-u) = ረ/'ሃ% - ም) +(e -&)  IV-10 
The iteration begins by substituting the approximations, so and Fy, 
into the right hand side of equation IV-10 and then invoking the condi- 
tions for periodicity. 


The most expedient way to do this is to write € as 


é-a 44a, GP. Ivi 


The value of / is then obtained from equation IV-3 and the result 
cubed. The value of m is then differentiated directly. This proce- 
dure eliminates the tedious progress of computing zn and then multi- 
plying by ያ : 

From equation IV-3, using equations IV-9a and IV-9b there results: 


ም -[£ p -ፉ 36 (178)] nin we 


T Bw La (45 Lio Sez. IV-12 
Letting 
p= # [ (4-8) 4 m (A 5፻)] IV-13a 
and 
a a መ B^), IV-13b 


there results upon cubing IV-l2: 


ፖ- HE + FG Janet -3 (56 +G Dee 
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3 | 
TF(E+ FG-L° )ansut 
+ ክክ ተ ባን yw č 

= # em PAT. IV-1l 


The condition for periodicity hence becomes 
3 2 

A-D w AF a ረ በር = ሪጋ, IV-15 
(A-B) ኘህ ህተ => w (4 e) 


The terms F and G are replaced by equations IV-13a and IV-13b 


respecitvely and upon simplification there results: 


Za lA-B)+ (17 Lg NAZE) + SM (AE) 


E lm ^(A-B + £3 Imm*(A1-8°)=0 


IV-16 
where 
AN == > a IV-17a 
ZZ z 2 IV-17b 
P = La IV-17c 
Se > ሙ / IV-17d 
lo =W FE | EC # IV-17e 
an = ፡1 = Æ IV-17f 
በቪ 

Ww? = - = wer IV-17g 


The result of substituting equations IV-9a and IV-9b into IV-6 is: 


Tz =4(AtB) + E5048?) was 












E ፣ 


- c — -] e» 








This is the equation for 4 found in chapter III. Omitting the 
trivial solution, A = B, equation IV-16 defines a solution curve, % ‘ 
Bistability results when E becomes tangent to ሪ . This occurs at a 
point on the A = B line. To find the value of Is required for bi- 
stability, equation IV-16 is first transformed into the polar coordin- 
ate form by the same transformations used in chapters II and III. The 
resulting equation is solved for Le , and this value is used to compute 
A. Equation IV-18 then produces the required Is. 


The transformations indicated above reduce equation IV-16 to 


4=0= -2 ፦ ታ(/#ፉ 2//#ሥ-#ዎያታ” 


7 
Il hrm Wl + ge mn flop S -— 


where ያማ and 4 (6) are defined as before, and 


L2 16 00 8-64 on ፉ 2ፖረሠሪ 22 ICO ፉ/..- ፲፻-20 

If the term, Ps is neglected, equation IV-19 is of the same form 
as the equation defining E . Since the functions involving the co- 
ordinate 4 , are always positive, the effect of the second non-linear- 
ity is to reduce the required current for bistability. The amount of 
reduction depends on the value of l,m, and n. 

There is no explicit statement derived, for the value which s must 
take. However, from previous developments of the single non-linearity 
circuit, it can be assumed that s5o, at least. 


Substituting the value O = O degrees into equation IV-19 results 
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“ጋና ‘ant 


Z £ I ያ 
This equation is a quartic ጨቃ? . The general form of the 


quartic is written as 


፲/ፆራ ደ x’ EZ CA += 2, IV-21 


where a,b,c and d are real, constant coefficients. 
To find the roots of equation IV-21, it is necessary to find any 


real root of the cubic: 


84 pe #8 9 ac to) = le * ሪ/2*- 24)] =, I-2 
The roots of equation IV-21 are then taken from the two quad- 


ratics: 


x [£-Yy2 »25-4]x *(q i eo 1-23 
L+ [E +y E g-t] alg (geo -2 


where 9 is the root found in equation IV-22. 

Obviously this is a tedious process. Fortunately, it needs to 
be done for only one value of O. However, to arrive at a general 
solution in terms of the parameters l,m,n, and s is not within the 
realm of this thesis. For any given components and source frequency, 
the value 6 «2 obtained through the quadric solution is substituted 


into the polar form of equation IV-11: 


5h 





420 = ES = +41 ($006 - (ea36)0 bateo ሪራ. 
Iv-2h 


This procedure gives the required value of Ea for bistability. 

In summary, the results of analysis of the doubly non-linear 
series-resonant trigger pair are as follows: 

1) The value of Is should be less for bistability in the double 
non-linearity circuit than for the single non-linearity circuit; 

2) There is no explicit development for the value which s should 
take, but it is implied by previous study that S 3 Ô, 

3) The required value of Ig is obtained from auxiliary quartic, 
cubic and quadratic equitem, as given by equations IV-20, IV-21 and 


IV-22. 
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CHAPTER V 
EXPERIMENTAL RESULTS 


The experimental results obtained with the circuits under consid- 
eration are given in this chapter. They are not offered especially to 
confirm theoretical developments, but to demonstrate that the circuits 
are operable. If the experiments agree with the theoretical develop- 
ments, so much the better. 

The charge-voltage relationships were determined through the use 
of the test equipment arrangement shown in figure V-1 for 60 cycle 
supply; and by the arrangement shown in figure V-2 for the 5 to 20 
kilocycle range. Sample hysteresigrams for frequencies of 60, 5000, 
10,000 and 20,000 cycles are shown in figure V-3. In figure V-l is 
the hysteresigram taken at 60 cycles, but with the vertical gain and 
intensity set so that the loop characteristics are distinguishable. 
The locus of the tips of a nest of 20,000 cycle major hysteresis 
loops is shown in figure V-5. 

For all experimental work the MUCCN type VSE non-linear capaci- 
tors were used. The values of the coefficients for these type capaci- 
tors were obtained from figure V-5 by approximating the average initial 
and final slopes. The Q-E relationship resulting from this procedure 
is shown as the dashed curve in figure V-5. 


The coefficient values thus obtained are 


a, >= #ያ/ሪ á volts/coul V-la 
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60 cycles C 


Figure V-1 





Test Equipment for obtaining hysteresigram at 
60 cycles. The VIVM should be peak reading. 





Figure V-2 


Test equipment for obtaining hysteresigrams from 


5 kes to 20 kes. The VTVM should be peak reading. 
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f: 20 kes | f: 10 kcs 
Emax: 350 volts Emax: 350 volts 
Qmax: 0.27 ucoul. Qmax: 0,3 ucoul, 





f: 5 kes f: 60 cycles 
Emax: 300 volts Emax: 375 volts 
Qmax: 0,299 ucoul. Qmax: 0.315 ucoul. 


(Q axis reversed) 


Figure V-3. Hysteresierems of MUCON type VSE non-linear capacitors, 





Figure V=-4. Hysteresigram of MUCON type VSE non-linerr capacitor 
showing open loop for a source frequency of 60 cycles. 
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V-1b 


V-le 


The circuit of figure I-la was first tested using the following 


component values: 










ES 
an 
20 kes 233 kes 233 ks 


For these components, the computed and experimental results are 


tabulated: 










s RE 31 1 2 1 45 M 
wi 
= 
La (exp) m 


Since the parameters, ay and a3, were not separately determined 





for the higher frequency, it is not suprising that the experimental 
and theoretical results disagree. However, from tests two and three, 
it was noted that the higher value of s gave a quicker "jump" but that 
the temperature of the crystal oven in which the capacitors were 
located rose 10 degrees centigrade suggesting a radical change of the 


capacitor parameters as a function of temperature. 
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Figure V-6, E, and E) for the two-branch series resonant circuit. 
The source voltage is larger for the waveforma on the right. 





Figure V-7, Ej and E5 for the two-mesh parallel resonant circuit, 
The source current, Iş, is larger for the waveforms on the right. 
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Figure V=8, Relative switching ratios and switching rates of the 
single non-linear resonant circuit and the double nonelinear cir- 
cuit, Figure a shows the switch for the two-branch series resonant 
circuit using a non-linear inductance; figure b shows the switch for 
the two-branch serier resonant circuit using both non-linear induct- 
ance and non-linesr oapacitance, 
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3000 pps 1000 pps 
300 pps 100 pps 


Figure V-9. Output voltage for various trigger retes of one branch 
of the two-branch series resonant circuit utilizing two non-linear 
elements. The voltage waveform is taken across one of the non-linear 


capacitors, The trigger pulse was 6.1 microseconds in duration and 
65 volts in amplitude, 








Switching was not obtainable by use of a pulse generator, but 
only by the use of a long-time (five secor’s or more) short circuit 
across the high current capacitor. This indicates that the capacitors 
were heating excessively and that the subsequent change in parameter 
values was detrimental to switching. Methods of maintaining both 
capacitors at the same constant temperature were not feasible because 
of limitation of time available. Figure V-6 shows the output voltages 
E and Eo at 233 kes after the bistability condition has been obtained. 
The phase relation conditions should be noted. The in-phase condition 
was first obtained, and by increasing Es, the out-of-phase condition 
resulted. 

The parallel circuit was next investigated at the frequency of 


20 kilocycles. The data on this test are given below. 









በው ውሉ 


Figure V-7 shows the output voltages, E, and E, after the bi- 






stability condition has been obtained. Note that the two voltages 
begin out-of-phase and approach the in-phase condition as Ig is 


increased. 


Losses in both the series resonant circuit and parallel resonant 
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circuit cause the large deviation from the 180 degree relationship. 
The fifth test concerned itself with .he double non-linear, two- 

branch series resonant circuit versus a single non-linear two-branch 

series resonant circuit. The nort-linear inductors were those used by 


Spitzer / 7 /. The parameters of these were: 


6, = 2 x 104 amps/oen V-2a 
ፌ - 2 x10! amps/ aen? V-2b 
ANS A Genus V-2c 
N = 200 V-2d 


The capacitor replacing the non-linear capacitance for the 
single non-linearity circuit was chosen to give approximately the 
same switching ratio, 


The results are given below: 
Nonlinear Capacitor 
Linear Inductor Nonlinear Inductor 


Figures V-8a and V-8b show the relative switching ratios and 














switching rates of the two circuits. The double non-linear circuit is 
shown to switch in about two cycles of the supply voltage whereas the 
circuit using only the non-linear inductance switches in about two and 
one-half cycles. 

The sixth test consisted of operation of the double non-linear 
two-branch series resonant circuit at various switching rates. The 


results of this are shown in figure V-9. A Hewlitt-Packard Model 212A 


6h 





Pulse generator was used as the trigger source. 

Only during a brief interval of time as switching obtained on 
the two-branch series circuit using a non-linear capacitor. It was 
apparent from the experiments that the temperature change of the capac- 
itors is the primary deterrent to switching. 

It is felt that the data presented confirms the theory of the 
first four chapters to a good degree. The design equations seem to be 
valid for operation near the frequency region in which the capacitor 
parameters were determined. At high frequencies, the losses are so 
great in comparison to the degree of non-linearity, that the circuit 


does not operate within the theory propounded by the analysis. 
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APPENDIX I 


TRIGGERING SCHEMES 


It was stated in chapters II and III that a circuit, in which 
controlled switching is obtainable between the two bistable points, 
would constitute a trigger pair. Obviously the circuit must also 
include means for obtaining this triggering. The purpose of this 
appendix is to present a few ideas as to how this may be obtained. 

The discussion will be limited to the suggested methods without too 
much thought as to their output requirements. One important question 
is left unanswered viz. by what method is memory produced. This is 
an important question, but is considered beyond the scope of this paper. 

When the pendulum was studied in chapter I, it was found that the 
division between periodic and aperiodic solutions was dependent on the 
energy of the system. The circuits studied in chapters I, II, and III 
were shown to pass from the unistable condition to the bistable condi- 
tion when a particular value of Es or Is was obtained. Although the 
implication was not forthright the approximation of the Q-E functions 
of these circuits to the sine function used in the pendulum case leads 
to the conclusion that the two states-unistable and bistable -are con- 
cerned with the amount of energy possessed by the individual branches 
or meshes. 


As an approximation, only the energy contributed by the funda- 
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mental component will be considered. 

Consider first, the two-branch series resonant circuit. There 
are two possible methods of trigger application: a pulse applied to 
the capacitor; or a pulse induced in the coil. In figure AI-1, both 
methods are shown. 


From equation II-37, the value of A for which bistability occurs 


A= at fis CE 
B 73 In ar sm AT-1 


and from equation II-39 the value of A for the maximum A in the bi- 


is given by 


stable condition is given by: 





The current in one branch is Wh and for each of the two states- 


unistable and bistable - it is given by: 
Z = co À I 


Ze = @& Amax AI-3b 
Since the energy is associated with the magnetic field for part 


of the cycle of the driving frequency, the energy contained by the 


branch can be written: 


EE #ረ/ = LA AT-ha 


for the unistable state, or 


2 
e dM LL 
for the bistable state. 
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DO ODO 


Trigger 


1 


X 


Trigger 
ፐ 


Figure AI-1 Two Trigger Schemes For Two Branch 
Series Resonant Circuit. 


H OO HR +3 





Figure AI-2 Two Trigger Schemes For Two-Mesh 
Parallel Resonant Circuit 
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The energy difference is thus given by: 


urb s ap 
Ue U = 3m (s ane ) AI-5a 


by using equations AI-l and AI-2. This is the amount of energy that 
this branch must lose before it is possible to switch it into the lower 
charge state. 

Since it was assumed that B ag ; the amount of energy this cir- 


cuit must gain by virtue of a pulse is given from equation AI-lja as 
D= 24 (s - gE) AI-6 
cm a, s " 
When branch "1" has lost the amount of energy specified by equa- 
tion AI-5a and branch "2" has simultaneously gained the amount speci- 
fied by equation AI-6a, switching can occur. In order to assure 


switching, the equality should be replaced with the unequality sign 


as follows 


Ly > wh (s - 24°) 


! AT-5b 
= w?L Lp 
|| Sau (S x: <=. ያም 


Since the trigger pulse is applied to both circuits, it must 
cause this exchange of energy to be as specified. At this point the 
investigation should inquire about memory. Although the subject 
memory is beyond the scope of this paper, a short attempt will be made 
to explain how the energy exchange may be caused. 

The trigger pulse adds energy to both branches, "1" and "2". The 


addition of energy to branch "1" does not change the effective capaci- 
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tance of the non-linear element by much since it was already operating 
in the saturated region. However, the energy addition to branch "2" 
does cause a large change in the effecting capacitance of the non-lin- 
ear element. Hence the opposition to the fundamental component de- 
creases to the point of non-linear resonance in branch "2" but does 
not decrease appreciably in branch "1". Since the source is provid- 
ing a constant current, the current now begins to divide between 
branches "1" and "2", adding to branch "2" by taking from branch "1". 
Branch "2" successively offers less opposition to the current whereas 
branch "1" successively increases its opposition. The process cascades 
until switching has occurred. When the energy in branch "1" has 
reached the value specified by equation AI-la, the trigger pulse should 
cease. Otherwise it would be possible to have branch "1" return to its 
original state of heavy conduction. The reason for this is that the 
pulse is supplying a certain fixed amount of energy which must divide 
between the two branches. Since the effective capacity of the unsat- 
urated branch is greater than that of the saturated branch, the unsat- 
urated branch is more receptive to the energy. However, when the 
branches reach the energy level specified by equating AI-la, the recep- 
tion is equally well for both branches. Only the effect of cascading 
can then cause the switch. 

The least amount of energy that the pulse must give is that speci- 
fied by equation AI-6a. Since the energy will divide, andat some time, 


equally, perhaps the pulse energy required should be: 


CE 





WR * 

= > Y S — ars ). 

Up aU መብ AI-6c 
The required energy expression for the two-mesh parallel reso- 


nant circuit can be similarly obtained. 
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APPENDIX II 


DERIVATION OF THE SECOND ITERATION SOLUTION CURVES 


The purpose of this appendix is to show how the second iteration 


solution curves were obtained. The second iteration begins with the 


approximations: 5 - 
fo - Aanwt DAT sut "m 
ም = B avnwt DE nA. AII-lb 


These are substituted into the right hand side of equation II-6, which 


after omission of the trivial solution, A=B, results in: 


Z0 = -5+ 041818?) 1m VHB AB YAB YE) 


Aon D AA B ABB ABEABIM 1-2 


First, the linear transformation: 
ለ] | -a|s 
" AII- 
B. | % 3 
is applied. This results in 
2 
P=0=- u +d (3K 78") + L(5 via ye) 
PEN Vioc) — AILh 
+2 (7 HE AB 8°): 
This equation is next transformed into polar coordinate form 


by the relations se Coad and Barmen 0 . This results in 


(see appendix III for expansion formulas): 
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E P 
bso- SrtA iso got ans 
where 
foe) = 2+ C698 ag 
46) = 3,5% 206026 - 0.5 080 46 AIT-6b 
Ale) = 5 + lsa 26/4- /200020 -Co2 20] 4- 799 u ee 


Equation AII-5 is now reduced to a cubic through the relation 


PER: This results in: 


#/)፡ሪ። -ሯ # EON + BIW IL HOT? AII-7 

From equations AII-6a, AII-6b, and AII-6c, it is noted that, 

for all 6, SC #)>0, a9) 20 689. #/#) > ጋ. . Thus the terms of 
equation AII-7 varying in sign only as does Jl. 


The coefficient signs are 
WA) a ት ጠታ 
tod Qo 
Thus there is one or no possible positive roots; and there are 
two or no possible negative roots. Since there are three roots, it 
can be stated that there are two negative roots and one positive root. 
The interest lies entirely with the positive root since the negative 
roots produce complex values off which have no meaning as far as the 
circuit is concerned. The "value! /l!-oeO does not satisfy the equation. 
Iherefore Ya is a single closed curve. 


It now becomes necessary to find the real root of Yin). 
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Ihe general form of a cubic equation is 
E ጃራ u +7 = 0 
D qns E Xn = ረ. ል፲፲-6 
The trigonometric method of solving cubics is used and the 


_ ህ(“ሙሑሙ 
AT AH E ) AII-9a 


real root is given by 





where 
ግዜ - 1455 Ca [5 ws in b 
241° 7 y-2* AII-9 
if ያሩ Bd 745^ $0] 
ረ... ፈይ ከክ. 
- iff €Ó and ያታ ጋ መመሪ) 
ዕ 
. i) 
24-215 2 JE 2 À] ል፲፲-94 
if f 20; 
and E AC — P AII-9e 


Z = ፊ /2ፈሪሪ - ad) -47 AII-9f 


The upper sign is chosen in the ፖ evaluation 4ያሪ4>ሪ and the 
lower sign if he O. It turns out that for all 0, 970 and thus 
equation AII-9a is the vehicle by which the roots of 7 are obtained. 


The equations of AN) are tabulated below for values of Q. 


ፀ Ws (r) 

0° -s/m 41.5 r 4125D rê 41.75 D*r2 
10° -8/8 4hr 41.00D x2 41.86 Der? 
20° -s/m 41.385r 41.23D r^ 42.065D^r2 
309 -s/m 41.25 r 41.19 Dr2 42.22 Der 
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9 Wr) 


1,00 -s/m 41.09 r 4 1.07 D r^, 2.15 Dr? 
15° -s/n 4 r 4100 Dr^ 2 D? 
509 -s/m 40.913r + 0.905 D r?4 2.07 D?r2 
60° -s/m 40.75 r 4 0.687 D r24 3.280ጭግ 
708 -s/m 40.617r 4 0.770 D ea 0.753D°r> 
80° -s/m 40.53 r 4 O.SID r?^ 4 0.37 Dr? 
900 -s/m 40.5 ፥ ፋ 0258 5: + 0.25 ፻ 


From these equations, the value of eis obtained for each angle 


specified 
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These are the only values needed to form a plot, since f4 is 
symmetrical about both the o< axis and 8 axis as can be seen from 
equation AII-L. 

The plots in figures II-l and III-2 were obtained from the 


above tabulations for an s value of 0.68. Although not labeled, 


the coordinates are in units of m ^. 


Ihe development of fg and p is done in a similar fashion. 
Tabulated equations of # ; for values of O are given below. Since A 
differs from # but little, except for large values of 2 or m, it 
has not been tabulated. The only real interest lies in the fact 
that both A and K have only one real positive root. For the value 
O - 90 degrees, all coefficients of ye vanish. Hence, # and £ are 
open curves. Since they were derived from monotonically decreasing 
functions, €, (9) and eala) by a linear substitution, they will 
also be monotonically decreasing functions about the axis of symmetry. 


Hence the shape is as shown in figures III-l and 111-2. 
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Using the method previously described for finding roots of the 


cubic, there arise solutions such as 


x a A orah /ሄመ” ረ መዴ), AII-10 
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where b and c are constants determined for each value of 9. For 


values of Calm Es <<] , the approximation 


ሪጋ talar a መሪ 27) ATI-11 


can be used. Since the approximation specifies x «| , there is the 


series equality: 


"Tu pU EU 
LH (X FI) = X-EK FER LX + = A II-12 


Using only the first term of AII-12, equation AII-10 becomes 
A ራ MEX ig fon Es ነ 
ሓ AN Y» 8 ZE AII-13 


Again since Cres Z<|, the hyperbolic sine function can 
be replaced by the first term of its power series expansion. Thus 
X, 7 ፌር AII-1) 
For the values of m and a used in experiments, Es should be 
less than 100 volts for all the approximations to be valid. It was 
on this basis that # and A were plotted in figures III-L and 
IIT-2 respectively. It is not anticipated that a marked difference 
would result if s value of Es were such that equation AII-10 must 
be used. 
The slope of the 4% curve is obtained by application of the 
formula 
LnT = Lama ARE f 
pa # መሬሪሠጋጋሪ ATI-15 
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where T is the angle the tangent makes with the horizontal axis. 
Since the horizontal axis is the œ axis, for the tangent to be 
parallel to the £ axis, the denominator of AII-15 must vanish, or 
the numerator must become infinite. 


Differentiating AII-5 with respect to 9, there results: 
ፊ ፖሪ። / J+ ታህ ap] y 
. | $2029 2% (Ojo wa AII-16 


From equations AII-6a, AII-Ób and AII-6c, it can be seen that 
£1(0), g'(0) and h'(0) each have the factor sin 9, sin 20, or sin 
ho. Thus, J vanishes for O - O degrees and therefore so does the 
denominator of AII-15. It has then been proven that the tangent of 


p at zi is parallel to the axis. 
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APPENDIX III 


MATHEMATICAL FORMULAS USED IN ANALYSIS 


This appendix lists, for the readers’ reference, mathematical 
formulas used in the analytical portion of this paper. 

The trigonometric expansions are obtained by applying the bi- 
nominal expansion to the exponential form of the sine and cosine func- 


tions. Thus there results 
- 3/m 
mx o 


ሠ m x) = ኤር. AIII-la 


am m nm 
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en am 


(uam x)= v ል፲፲፲-25 
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ን 2m+l / 
kunna) = - ያ” <er (Qml)m X 
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In determining that Y was an ellipse and investigating if it 
became an imaginary ellipse due to the presence of losses, the follow- 
ing properties of conic sections were used. 


The general form of a conic is written as 


4 2" ፉሪ2#2 ሥ፡ያ ሥ22። “ዖ/ሪ ፇር. uns 


The relations, D, d and c are defined as: 


a ከ ጀ 

D = h b f AIII-6a 
g f c 
a h 

d = ATIT-6b 
h b 
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e ee a+ b ALII-6c 

The cónicas: 

l. a parabola if d - O, D 0; 

2. parallel lines (may be coincident or imaginary) if 

d=0, D=0; 

3. an ellipse if d>0, eD<O; 

lh. no locus (imaginary ellipse) if d >0, eD>0; 

5. a point ellipse if d>0, D=0; 

6. an hyperbola if d<0, D¥O 

7. two intersecting lines if d« O0, D= 0. 
In order to remove the xy term of the generalized conic equation, 

the rotation matrix 


የዳ መማ) 


ሮ| 1ጩሬሪ  ርሪሪ!| |ዶ — 


is applied, where 9 is given by 
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